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Abstract 

Beginning  with  fundamental  properties  such  as  passivity  or 
incremental  passivity  of  the  network  elements  comprising  a 
switched  power  converter,  the  nominal  open- loop  operation 
of  a  broad  class  of  such  converters  is  shown  to  be  stable  in 
the  large  via  a  Lyapunov  argument.  The  obtained  Lyapunov 
function  is  then  shown  to  be  useful  for  designing  globally 
stabilizing  controls  that  include  adaptive  schemes  for  han¬ 
dling  uncertain  nominal  parameters.  Numerical  simulations 
illustrate  the  application  of  this  control  approach  in  DC-DC 
converters,  v 

1  Introduction 


[12]  and  [14]  also  consider  control  design  using  Lyapunov  func¬ 
tions. 

The  paper  is  organized  as  follows.  Section  2  develops  fun¬ 
damental  stability  properties  for  a  broad  class  of  switching 
converters.  In  particular,  we  derive  the  form  of  a  Lyapunov 
function  that  illustrates  that  each  member  of  this  class  is 
open- loop  stable.  In  Section  3,  the  Lyapunov- based  control 
approach  is  introduced.  This  section  includes  a  simple  ex¬ 
ample  to  demonstrate  the  method,  an  outline  of  the  general 
approach,  and  an  illustration  of  how  an  adaptation  scheme 
can  be  incorporated  to  handle  uncertainties  in  the  nomiual 
operating  point. 

2  Open-Loop  Stability  of  Switch¬ 
ing  Converters 


Most  control  schemes  for  power  electronic  circuits  in  present 
use  are  obtained  by  linearizing  a  nonlinear  model  about  a 
nominal  operating  point  or  trajectory.  Large  signal  transients 
that  occur  at  power  up  or  overload  recovery  are  handled  in 
an  ad  hoc  manner.  Designers  analyze  each  circuit  individ¬ 
ually  to  prescribe  a  scheme  to  accomodate  a  designated  set 
of  large  signal  transients.  This  paper  addresses  the  issue  of 
how  one  might  do  better.  In  particular,  the  paper  develops 
a  methodology  for  designing  control  laws  for  fast-switching 
converters  that  result  in  globally  stable  behavior,  robustness 
against  parametric  uncertainty,  and  satisfactory  transient  re¬ 
sponse.  The  main  approach  in  this  paper  is  based  on  the  use 
of  Lyapunov  functions. 

There  has  been  considerable  previous  work  along  these 
lines.  The  paper  of  Erickson  et  al.  [11]  develops  a  large-signal 
averaged  model  for  switched  converters  and  points  out  the 
hazards  of  control  designs  based  upon  small-signal,  linearized 
models.  The  previous  work  on  huge  signal  control  schemes 
can  be  divided  roughly  into  two  groups.  One  is  based  on  so- 
called  “switchiug-law"  controls  where  the  position  of  a  con¬ 
trolled  switch  is  directly  commanded  as  a  function  of  the  in¬ 
stantaneous  values  of  the  circuit  variables.  Examples  of  these 
are  the  sliding  mode  control  schemes  of  [1,  2,  3,  4,  5,  6,  7,  21] 
(which  include  current-mode  control)  and  the  bang-bang  con¬ 
trol  schemes  of  [9, 10, 15].  The  second  approach  relies  on  the 
state-space  averaged  model  for  the  converter  of  interest.  Rep¬ 
resentative  schemes  are  those  in  [12,  13,  14,  22],  References 
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In  this  section,  switching  converter  systems  (understood  to 
include  source  and  load)  that  consist  of  an  interconnection 
of  ideal  DC  sources,  ideal  switches,  incrementally  passive  re¬ 
sistors,  and  passive  linear  reactive  elements  are  considered. 
Multiport  circuit  elements  are  included  iu  the  development 
here.  This  class  of  switching  converters  is  shown  to  l>e  stable 
by  exhibiting  a  Lyapunov  function  that  corresponds  to  the 
energy  in  the  increment  with  respect  to  au  arbitrary,  nominal 
state  trajectory.  The  argument  is  exteuded  to  include  non¬ 
linear  reactive  elements  that  are  strictly  relatively  passive  in 
the  case  where  the  switching  frequency  becomes  infinite,  and 
stability  with  respect  to  an  equilibrium  point  is  considered. 
Essential  background  on  network  theoretic  issues  for  the  de¬ 
velopment  in  this  section  is  contained  in  Appendix  A.  (Also, 
see  [16, 17,  18].) 


2.1  Switching  Converter  Stability  Under 
Finite  Switching  Frequency 

Let  the  switching  converter  be  composed  of  ideal  DC  sources, 
ideal  switches,  incrementally  passive  resistors,  and  linear  pas¬ 
sive  reactive  elements.  A  diode  may  be  considered  as  either 
an  ideal  switch  or  as  an  incrementally  passive  resistor.  We  Q 
suppose  the  switches  are  operated  in  accord  with  a  given  ar-  Q 
bitr&ry  switching  pattern  and  suppose  that  we  arc  given  a 
nominal  solution  corresponding  to  the  given  switching  pat-  — 
tern. 


For  each  branch  of  the  network,  denote  the  nominal  tra-  _ 

jectory  by  (f>((),  «(<)},  and  form  the  (not  necessarily  small) 
increments  with  respect  to  the  nominal  trajectory  for  each  _ 

Jes 
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network  branch,  i.e. 


6i(t)  =  i(t)  —  i(t) 

£v(t)  =  v(t)  —  v(t)  (1) 

By  applying  Tellegen’s  theorem  to  the  increments  in  all  the 
network  branches  when  the  circuit  is  in  any  one  of  its  topolo¬ 
gies,  we  obtain 

0  =  52  H6v+  52  + 

DC  »o*tr  c cj  Sw  tick*  i  Res. 

Y^6i6v+'E6i6v  (2) 

Ind.  Cep. 

The  summation  involving  DC  sources  is  always  zero  since 
the  increment  in  either  voltage  or  current  of  each  term  is 
necessarily  zero.  The  terms  involving  switches  also  add  zero 
contribution  to  the  sum  in  (2)  for  the  same  reason.  The  third 
summation  on  the  right-hand  side  of  (2)  is  always  nonnegative 
since  each  term  is  individually  positive  or  zero  as  a  result 
of  the  incremental  passivity  of  the  resistors.  In  fact,  this 
summation  can  be  thought  of  as  the  dissipated  power  in  the 
increment  with  respect  to  the  given  nominal  trajectory.  The 
remaining  two  summations  represent  the  time  rate  of  change 
of  the  stored  energy  in  the  increment  Combining  these  facts, 
we  can  write: 

^-V(6x)  =  52  &*&v+  ^siSv^-Y,  (3) 

/  nd.  Cap.  Ret. 

where 

V(<X)  =  52(l/2)(«i*)*I4(*i*)  + 

led. 

DO/2 ■)(**)•<:*(*»).  (4) 

Cap. 

(Note  that  the  superscript  *  indicates  transpose.)  Because 
of  the  assumption  on  passivity  of  the  reactive  elements,  the 
quantity  V(6.r)  which  we  shall  from  now  on  refer  to  by  the 
suggestive  name  energy  in  the  increment,  is  a  positive  definite 
quadratic  function  of  the  incremental  state  variables.  Since 
equatiou  (2)  holds  identically  for  any  of  the  possible  circuit 
topologies,  the  energy  in  the  increment  is  a  Lyapunov  func¬ 
tion  for  the  dynamical  system.  In  particular,  we  have 

-^V(Jx)  =  -  52  <  0-  (5) 

In  conclusion,  the  energy  in  the  increment  is  a  Lyapunov 
function  for  the  given  nominal  trajectory,  and  we  see  that  the 
nominal  trajectory  is  stable  in  the  large.  Since  the  nominal 
trajectory  selected  above  can  be  taken  as  any  solution  trajec¬ 
tory,  this  statement  implies  that  any  two  solution  trajectories 
do  not  diverge. 

Typically,  asymptotic  stability  in  the  large  can  be  con¬ 
cluded  as  well  since  at  least  some  parasitic  loss  is  always  asso¬ 
ciated  with  each  energy  storage  element,  i.e.  series  resistance 
with  inductors  and  parallel  leakage  resistance  with  capaci¬ 
tors.  An  argument  for  asymptotic  stability  appears  in  [16] 
for  circuits  that  have  a  DC  equilibrium  point,  and  consist  of 
only  two-terminal  elements.  In  [16],  lossiness  is  guaranteed  to 
be  associated  with  each  state  variable  by  excluding  inductor- 
capacitor- voltage  source  loops  and  inductor-capacitor-current 


source  cutsets.  (In  this  case,  we  would  require  that  all  resis¬ 
tors  be  strictly  incrementally  passive.) 

A  special  case  of  the  above  result  is  when  the  switches  are 
operated  with  a  periodic  switching  pattern,  and  there  exists  a 
nominal  periodic  steady  state  solution.  In  this  case,  the  result 
states  that  the  given  periodic  steady  state  trajectory  is  stable 
in  the  large.  This  result  is  of  particular  interest  for  the  case 
of  a  DC-DC  converter  operating  with  constant  switching  fre¬ 
quency.  Note  that  this  result  holds  up  for  DC-DC  converters 
operating  in  the  discontinuoas  conduction  mode.  This  can  be 
seen  by  redrawing  the  schematic  for  the  DC-DC  converter  of 
interest  with  an  ideal  SPDT  switch  and  incrementally  passive 
resistive  device  (i.e.  diode)  replacing  each  trims  is  tor -diode 
pair.  For  example,  we  would  redraw  the  up-down  converter 
of  Figure  2  as  shown  in  Figure  1.  The  circuit  of  Figure  1  sat¬ 
isfies  the  conditions  for  its  nominal  periodic  trajectory  to  be 
stable  in  the  large,  and  it  makes  no  difference  whether  or  not 
the  nominal  trajectory  contains  a  portion  where  the  inductor 
current  is  identically  zero. 
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Figure  1:  Up-Down  Converter  Redrawn  to  Dlustrate  Stability 
in  Case  of  Discontinuous  Conducion 


2.2  Stability  under  Infinite  Switching  Fre¬ 
quency  and  Constant  Duty  Ratio 

The  result  given  above  can  be  extended  to  the  case  where 
the  switching  frequency  becomes  infinite.  Infinite  switching 
frequency  actually  corresponds  to  the  state-space  averaged 
model  for  a  converter,  and  in  this  way  an  open-loop  stability 
result  can  be  obtained  for  state-space  averaged  models.  The 
main  difference  from  the  case  of  finite  switching  frequency  is 
that  one  needs  to  consider  the  stability  of  an  equilibrium  point 
for  an  averaged  model  of  a  DC-DC  converter  rather  than  a 
limit  cycle.  In  such  a  set-up,  it  is  possible  to  include  nonlin¬ 
ear  reactive  elements  as  well  as  nonlinear  resistive  elements. 
This  is  of  interest  in  power  electronic  circuits  since  nonlinear 
reactive  elements  do  occur  in  practice.  The  following  theorem 
summarizes  the  result  for  infinite  switching  frequency. 

Theorem  2.1  Suppose  that  o  switching  converter  is  con¬ 
structed  from  ideal  switches,  ideal  DC  sources,  incrementally 
passive  resistors,  reactive  elements  that  are  strictly  relatively 
passive,  and  that  its  averaged  model  has  an  equilibrium  point, 
then  the  equilibrium  is  stable  in  the  large. 

This  theorem  is  easily  proven  by  demonstrating  the  exis¬ 
tence  of  an  appropriate  averaged  circuit  model.  See  [23]  for 
more  details.  We  shall  rely  heavily  on  this  result  in  the  sequel 
where  control  laws  based  on  state- space  averaged  models  are 
developed. 
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3  Lyapunov-Based  Control  Design 

In  this  section,  an  approach  to  control  of  switching  power 
converters  based  on  the  use  of  Lyapunov  functions  will  be 
introduced.  The  main  focus  will  be  on  control  d"*sign  based 
on  the  state-space  averaged  model  for  a  given  switching  con¬ 
verter.  The  converters  of  interest  are  those  that  satisfy  the 
conditions  guaranteeing  that  nominal  state  trajectories  are 
globally  stable  under  open-loop  operation,  specifically  con¬ 
verters  constructed  from  incrementally  passive  resistors,  ideal 
sources,  ideal  switches,  and  passive  linear  reactive  elements. 
One  particular  choice  of  Lyapunov  function  for  centred  de¬ 
sign  purposes  that  will  be  of  interest  is  the  energy  in  the 
increment. 

We  shall  begin  by  illustrating  the  Lyapunov-based  control 
method  with  an  application  to  an  up-down  converter.  Then, 
we  shall  demonstrate  how  such  a  control  design  can  be  ob¬ 
tained  in  a  more  general  way.  There  is  typically  some  freedom 
in  the  choice  of  Lyapunov  function  for  the  control  design,  but 
we  shall  exhibit  some  particular  advantages  of  using  the  en¬ 
ergy  in  the  increment  Finally,  we  shew  how  &n  adaptation 
scheme  can  be  incorporated  to  handle  parametric  uncertainty. 
Generalizations  to  converters  containing  nonlinear  circuit  el¬ 
ements.  to  converters  that  handle  time-varying  input-output 
waveforms,  and  to  converters  operating  in  the  discontinuous 
conduction  mode  are  given  in  the  thesis  [23].  A  method  (dual 
to  the  control  design  approach)  for  designing  state  observers 
is  also  considered  in  [23]. 


7,  =  2amps 
om  =  —  9volts 
>a  =  3.2amps 

The  relevant  matrices  of  the  system  are  as  follows: 


A 

B 

b 

Q 


f  0  (i -*.)/!  1 

[-(l-da)/C  0  | 

0  -l/L 

1/C  0 

(V.  -  Va  )/L 
in/C 

fiSl 


The  result  on  open-loop  stability  in  Section  2  guarantees  that 
the  energy  in  the  increment  is  a  Lyapunov  function  for  open- 
loop  operation  of  this  circuit.  For  the  up-down  converter,  the 
energy  in  the  increment  takes  the  form 

V  =  iL(,-,.)>  +  ic(«.-r.)1,  (7) 

or 

V  =  l-r‘Qx.  (8) 

Differentiating  V  along  the  system  trajectories,  we  obtain 


3.1  Example:  Up-Down  Converter 


Figure  2:  Up-Down  Converter 

Consider  the  up-down  converter  of  Figure  2  which  has  a  state- 
space  averaged  model  of  the  form 

x'  =  Ax  +  (Bx  +  b)d,  (6) 

where  the  two-component  state  x  consists  of  the  deviation 
of  the  inductor  current  from  its  nominal  value  (x\  =  t  — 1») 
and  the  deviation  of  the  capacitor  voltage  from  its  nominal 
(xj  =  r  —  t>„),  and  where  the  input  d  is  the  deviation  in 
the  duty  ratio  from  its  nominal  value  (d  =  d,  —  d„).  (Note 
that  d,  indicates  the  total  duty  ratio  here.)  The  parameter 
values  listed  below  were  selected  for  operation  at  a  switching 
frequency  of  SOKHz. 

C  =  5.4/iF 
L  =  O.lSmH 
R  =  o o 
dn  =  3/8 
V,  =  15  volts 


jV(«)  =  ir-(Q.4  +  A-Q)x  + 

i{x*(QB  +  B-Q)x  +  2  b'Qx)d.  (9) 

It  turns  out  that  QA  +  A‘Q  —  0  for  this  example,  which 
verifies  that  the  energy  in  the  increment  is  a  Lyapunov  func¬ 
tion  for  open-loop  operation  (d  =  0).  In  this  example,  it  is 
also  true  that  QB  +  B’Q  =  0.  These  relationships  hold  be¬ 
cause  of  the  lossless  nature  of  the  example  converter,  i.e.  the 
lack  of  resistive  elements  in  the  converter.  Considering  these 
relationships,  (9)  simplifies  considerably  to 

iv(x)  =  (b'Qx)d.  (10) 

Mw  stabilizing  control  schemes  can  be  obtained  by  in- 
spcctK/u  of  (10).  We  shall  consider  the  simple  control  law 
d  —  V  7*  with  a  real  and  positive,  modified  to  handle  the 
dut>  3  saturation  constraint  —  dn  <  d  <  1  —  dH,  i.e. 

f  -<*y.  -dn  <  d  <  1  -  dn 

d  =  <  -dn,  -ay  <  -dm  (11) 

l  I  -  dn,  -cry  >  1  -  dn 

where  y  =  b’Qx.  Here,  the  variable  y  takes  the  form 

y  =  (V,  -  t>„)(«  -  i„)  +  iH(v  -  vn) 

=  (*'.  “  i’)(«  -  in)  +  »(v  -  «’’»)•  (12) 

Note  that  the  only  dependence  on  circuit  parameters  is  on  the 
nominal  values  of  the  inductor  current,  the  capacitor  voltage, 
and  the  source  voltage.  This  property  is  shared  by  analagous 
control  schemes  based  on  the  energy  in  the  increment  for 
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many  other  switching  converters,  u  will  be  discussed  in  Sub¬ 
section  3.3.  The  dependence  on  nominal  values  of  circuit 
variables  is  of  crucial  importance,  and  this  issue  is  addressed 
in  Subsection  3.4.  There,  a  method  for  adaptively  estimating 
these  values  is  developed. 

To  investigate  the  closed-loop  behavior,  we  examine  the 
derivative  of  the  Lyapunov  function  V(x)  along  the  closed- 
loop  system  trajectories: 

±V(x)  =  yd 

-ay*,  -d„  <  d  <  1  -d, 

-dnV,  -ay  <  -dm  (13) 

(1  -  dn)y,  -ay  >  1  -  dn 

In  the  saturated  regions  (the  second  and  third  lines  of  (13)), 
the  time  derivative  of  V(x)  is  strictly  negative  since  either 
V'  <  — odj  or  V'  <  —  a(l  —  dn)i.  As  a  result,  state  trajecto¬ 
ries  quickly  enter  the  unsaturated  region.  In  the  unsaturated 
region  (the  first  line  in  (13)),  V(x)  is  strictly  decreasing  if 
y  /  0.  and  asymptotic  stability  can  be  concluded  by  LaSalle's 
theorem  since  y  =  0  is  not  a  system  trajectory  unless  x  =  0. 
To  see  this,  note  that  y  s  0  implies  d  =  0  and  the  following: 

b'Qx  =  0 

b'QAx  =  0,  (14) 

with  the  last  line  in  (14)  obtained  by  noting  that  y’  —  0.  The 
existence  of  a  nonzero  solution  x  to  (14)  is  equivalent  to  the 
statement  that  the  pair  (b'Q,  ,4}  is  unobservable.  However, 
this  pair  is  observable  in  this  example,  and  therefore  there 
are  no  system  trajectories  that  do  not  converge  to  the  origin. 

In  itis  example,  we  have  not  considered  the  effect  of  lossi- 
ness  due  to  parasitic  and/or  load  resistances.  The  effect  of 
such  passive  resistances  would  only  enhance  our  stability  re- 
suit,  by  causiug  additional  nonpositive  terms  of  the  form 
-x’Rx  (with  R  positive  semi-definite)  to  be  added  to  the 
terms  on  the  right-hand  sides  of  (13). 


Figure  3:  Root-Locus  for  Linearized  Closed-Loop  Control 
System 

Asymptotically,  the  decay  of  the  Lyapuuov  function  V(x) 
is  con  Trolled  by  the  eigenvalues  of  the  small  signal  model  ob¬ 
tained  by  linearizing  the  closed-loop  system  about  x  =  0.  In 
this  example,  there  is  some  freedom  in  placing  the  eigenvalues 


of  the  linearized  closed- loop  system  by  choice  of  the  gain  o. 
A  root  locus  of  the  closed-loop  eigenvalues  of  the  small- signal 
model  is  shown  in  Figure  3.  To  minimize  the  maximum  of 
the  real  parts  of  the  eigenvalues,  for  example,  the  gain  should 
be  selected  so  that  the  two  eigenvalues  coincide  on  the  real 
axis  at  — 20.05Krad/sec.  An  easy  calculation  indicates  that 
the  value  of  the  gain  required  to  obtain  this  eigenvalue  place¬ 
ment  is  approximately  a  =  .00785.  In  the  remainder  of  the 
discussion  of  this  example,  a  value  of  the  gain  of  a  »  .008 
will  be  used.  The  resulting  closed-loop  eigenvalues  are  at 
about  —  16.7Krad/eec  and  —  24Krad/sec.  Note  that  in  this 
example  the  dynamical  behavior  of  the  small  signal  closed- 
loop  dynamics  is  limited  by  the  natural  resonant  frequency 
(1  -  dn)/VLC  of  the  open-loop  state-space  averaged  system. 
Since  the  bandwidth  of  the  closed-loop  dynamics  is  usually 
designed  to  be  approximately  an  order  of  magnitude  below 
the  switching  frequency,  and  since  this  is  also  a  typical  res¬ 
onant  frequency  of  the  open-loop  dynamics  for  a  reasonably 
designed  converter,  the  preceding  limit  on  attainable  closed- 
loop  bandwidth  is  acceptable. 

We  expect  the  closed-loop  system  to  be  very  well  behaved, 
and  this  is  confirmed  by  the  digital  computer  simulation 
shown  in  Figure  4.  In  the  following  subsection,  we  present 
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Figure  4:  Digital  Simulation  of  Up-Down  Converter  under 
Nonlinear  Feedback  Control  Scheme 


a  derivation  of  a  class  of  control  schemes  to  which  the  alx>ve 
example  belongs. 

3.2  A  Basic  Approach  to  Lyapunov-Based 
Control  Design 

In  this  subsection,  we  show  how  to  derive  a  class  of  centred 
laws  for  a  switching  converter  model  of  the  form  (6),  to  which 
the  example  (11)  belongs.  Note  that  the  open- loop  stability 
of  the  system  (6)  is  crucial  for  this  approach,  and  hence  we 
shall  restrict  attention  to  switching  converters  that  satisfy  the 
conditions  guaranteeing  stability  under  nominal  duty  ratio 
operation.  A  basic  first  step  in  this  approach,  as  illustrated 
above,  is  the  specification  of  a  Lyapunov  function  for  open- 
loop  operation.  The  model  (6)  is  linear  and  time-invariant 
in  the  case  of  open-loop  operation  under  a  constant  nominal 
duty  ratio,  i.e.  d  =  0.  Since  the  open- loop  model  is  known 
a  priori  to  be  stable,  it  is  generally  possible  to  determine 
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a  family  of  suitable  quadratic  Lyapunov  functions.  In  fact, 
in  the  case  where  the  matrix  A  is  asymptotically  stable,  it  is 
possible  to  parametrize  the  family  of  such  quadratic  functions 
with  tbe  Lyamiuov  equation 

A-Ql  +  QiA  =  ~P,Pl-  (15) 


where  {/*,*,  A)  is  an  observable  pair.  The  existence  of  a  pos¬ 
itive  definite,  symmetric  solution  (Ji  to  (15)  is  guaranteed 
by  the  stability  of  the  matrix  A  and  the  observability  of  the 
pair  {Pj’M}  [8]-  See  [23]  for  a  method  of  selecting  a  suit¬ 
able  matrix  Q i  for  the  case  where  the  matrix  A  has  (simple) 
eigenvalues  on  the  jui-axis. 

Having  determined  the  form  of  a  suitable  matrix  Q  j,  it  is 
straightforward  to  specify  a  globally  stabilizing  control  law  for 
the  model  (6)  of  the  form  (11),  but  based  on  the  Lyapunov 
function  V'(x)  =  %x’Q\X,  as  follows: 


d  = 


-ay,  -dn  <  -oy  <  1  -  dH 
-dn,  -ay  <  -dH 
1  -  dn,  -oy  >  1  -  dn 


(16) 


where  y  =  (Bx  +  5)*Cjix.  One  particular  choice  for  Q\  is  Q 
where  I'(x)  =  i.r"Qr  is  the  energy  in  the  increment,  and  it 
turns  out  that  tnis  choice  leads  to  certain  nice  features,  which 
are  elaborated  below. 


3.3  Advantages  of  the  Use  of  the  Energy 
in  the  Increment  for  Control  Purposes 

As  noted  in  the  previous  subsection,  there  is  typically  some 
freedom  in  the  choice  of  the  Lyapunov  function  that  can  be 
used  in  the  control  designs  described  there.  Here,  we  outline 
three  advantages  obtained  by  using  the  energy  in  the  incre¬ 
ment  as  the  Lyapunov  function  in  these  control  schemes. 

Our  advantage  of  the  choice  of  the  energy  in  the  increment 
as  the  Lyapunov  function  for  control  design  purposes  arises 
in  the  computation  of  the  variable  y  =  (Bx  +  b)’Qx  which  is 
used  in  these  control  schemes.  In  particular,  one  can  always 
(indirectly)  measure  the  vector  Q(Bx  +  6).  To  see  this,  con¬ 
sider  the  modification  of  (6)  where  we  multiply  this  equation 
on  the  left  by  the  matrix  Q,  giving 


Figure  5:  Up- Down  Converter  with  Additional  Filter  Sections 


performing  the  described  measurement  process,  it  is  possible 
to  obtain  an  accurate  measurement  of  the  vector  Q(Bx  +  6). 
Consequently,  one  can  compute  the  variable  y  =  x’Q(Bx  +  b) 
by  forming  the  inner  product  of  x  and  Q(Bx  +  b).  The  only 
parametric  dependence  is  therefore  on  the  nominal  state  val¬ 
ues  required  to  determine  x,  the  deviation  in  the  states  from 
their  nominal  values. 

In  certain  cases,  it  is  possible  to  further  runplify  the  mea¬ 
surement  of  Q(Bx  +  b).  For  these  cases,  it  is  possible  to 
directly  measure  the  vector  Q(Bx  +  b )  by  measuring  certain 
branch  voltages  and  branch  currents  in  the  circuit  at  one  time 
inst  ant.  One  such  example  is  the  up-down  converter  of  Figure 
2.  For  this  example,  the  vector 


Q(Bx  +  b)=  V‘." 


(18) 


and  these  quantities  can  be  directly  measured  on  the  con¬ 
verter  circuit.  Necessary  and  sufficient  conditions  for  such  a 
simple  measurement  of  the  vector  Q(Bx  +  b)  are  given  in  [23]. 

A  second  potential  advantage  of  the  choice  Q|  —  Q  in 
(16)  is  that  it  is  possible  to  use  a  nearly  linear  version  of 
this  control  algorithm  by  replacing  y  =  (Bx  +  b)'Qx  in  (16) 
with  yt,„  =  b‘Qx,  and  still  maintain  global  stability.  (Of 
course,  the  saturation  constraints  are  still  in  effect.)  To  see 
that  global  stability  is  maintained,  consider  the  following  Lya¬ 
punov  analysis  with  V(x )  =  |x'Qx: 


jV(x)  =  ix*((.4  +  dB)'Q  +  Q(A  +  dB)]x  +  (b'Qx)d.  (19) 


Qx'  =  QAx  +  Q(Bx  +  b)d.  (17) 

Now  tbe  vector  on  the  left-hand  side  of  (17)  is  composed  of  the 
time  derivatives  of  the  inductor  fluxes  and  the  time  deriva¬ 
tives  of  the  capacitor  charges.  The  elements  of  this  vector 
are  necessarily  inductor  voltages  and  capacitor  currents.  The 
vector  Q(Bx  +  b)  is  the  amount  by  which  this  vector  changes 
when  the  duty  ratio  steps  from  —dn  to  1  -  rf„ ,  or  equivalently, 
the  amount  this  vector  changes  when  the  switch  configuration 
is  changed.  In  general,  it  is  possible  and  feasible  to  determine 
the  vector  Q(Bx  +  6)  during  each  cycle.  To  do  this,  for  each 
inductor  branch  one  would  measure  the  voltage  across  the 
branch  in  each  of  the  two  switch  configurations,  and  then 
form  the  difference  of  the  two  measurements.  This  difference 
constitutes  the  element  of  Q(Bx+b)  corresponding  to  the  par¬ 
ticular  inductor  port.  In  the  case  of  a  capacitor,  one  would 
measure  the  current  flowiug  into  the  capacitor  in  each  of  the 
two  switch  configurations,  and  form  the  difference  of  the  two 
measured  curreuts.  This  difference  constitutes  the  element 
of  Q(  3x  +  l)  corresponding  to  the  particular  capacitor.  By 


Now  the  first  term  on  the  right-hand  side  of  (19)  is  always 
nonpositive.  This  follows  from  the  fact  that  the  energy  in 
the  increment  takes  the  form  \x‘Qx  for  any  nominal  duty 
ratio,  with  the  fixed  matrix  Q.  The  choice  of  the  control  in 
(16)  (using  yim)  forces  the  second  term  on  the  right-hand  side 
of  (19)  to  be  nonpositive.  Global  stability  results  from  the 
nonpositivity  of  the  right-hand  side  of  (19).  Hence,  the  choice 
of  C?i  =  <?  in  (16)  permits  the  use  of  a  feedback  control  that 
requires  only  tbe  computation  of  the  linear  variable  ylin. 

A  third  advantage  of  the  use  of  the  energy  in  the  increment 
as  a  Lyapunov  function  for  control  design  is  that  a  control  law 
of  the  form  (16)  with  Qi  =  Q  can  result  in  global  stability  of 
a  more  complex  power  system  in  which  the  original  converter 
is  embedded.  In  particular,  if  the  converter  is  interconnected 
only  with  (relatively)  passive  circuit  elements,  the  resulting 
interconnected  system  is  always  guaranteed  to  be  stable.  For 
example,  if  an  additional  section  of  output  filter  is  added  to 
the  up-down  converter  of  Example  2,  as  shown  iu  Figure  5, 
the  control  law  designed  for  the  original  converter  stabilizes 
tbe  modified  circuit. 
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3.4  Adaptive  Control  Method  to  Handle 
Uncertain 

Nominal  State  Values 

This  subsection  considers  a  control  design  of  the  form  (1 1)  for 
the  model  (6),  but  in  the  case  where  the  nominal  state  vector 
is  unknown.  The  effect  of  this  uncertainty  is  to  replace  the 
variable  y  by 

y  =  (Bx  +  b)'Q(x  -  Sxm) 

where  (xm  is  the  uncertainty  in  the  nominal  operating  point, 
that  is 

Sxn  =  i,  -  r.  (20) 

where  r.  is  an  estimate  of  the  nominal  operating  point. 

To  implement  the  self- tuning  scheme,  we  shall  include  as 
part  of  an  augmented  state  vector,  an  estimate  x„(t)  of  the 
constant  nominal  value  of  the  state  vector  for  the  original 
plant.  We  can  equivalently  represent  this  estimate  by  its  er¬ 
ror,  i.e.  =  i„(f)  —  The  update  law  for  Sx„(t)  is 

selected  by  considering  the  Lyapunov  function 

V  =  ±x'Qx+±(6xn)-K(6xn)  (21) 

where  K  is  a  symmetric  positive  definite  matrix  and  Q  is  as 
previously  specified.  In  particular,  it  is  possible  to  stabilize 
the  system  by  choosing  the  update  law  to  be 

jt(6xn)  =  -K-'Q(Bx  +  b)  (22) 

in  conjunction  with  the  control  law  (11).  Note  that  y  can 
now  l>e  determined  without  any  uncertainty  arising  from  the 
unknown  nominal  state  values  since 

x-Sxn  =  (i,  -  *„)  -  (f„  -  *„)  *  x,  -  i„, 

where  x,  is  the  actual  full  state  value  (which  can  be  measured) 
while  is  stored  in  the  controller.  Note  that  it  is  generally 
possible  to  obtain  an  accurate  measure  of  Q{Bx  +  6)  as  dis¬ 
cussed  in  Subsection  3.3. 


Note  that  in  this  model,  the  quantities  without  subscripts  are 
deviations  from  nominal,  the  quantities  with  subscript  t  are 
total  variables  that  can  be  measured,  and  the  quantities  with 
subscript  n  are  nominal  variables.  We  only  attempt  to  esti¬ 
mate  the  nominal  inductor  current  since  the  other  nominal 
state  variable  (the  capacitor  voltage)  is  known.  The  output  y 
of  this  model  can  be  determined  exactly  since  i  -  (Sin)  is  pre¬ 
cisely  it  —  »«,  i.e  the  difference  between  the  actual  inductor 
current  and  the  present  estimate  of  the  nominal  value  of  this 
current.  The  control  design  can  be  completed  by  specify¬ 
ing  k  >  0  and  a  feedback  gain  a.  These  parameters  may 
be  selected  by  considering  the  small  signal  behavior.  For 
example,  with  a  nominal  load  current  of  2amps,  the  eigen¬ 
values  of  the  small  signal  linearized  model  ran  be  placed 
at  —7.713  ±  jl2.9Krad/sec  and  -11.36Krad/sec  by  select¬ 
ing  k  =  2778  and  the  unsaturated  gain  a  —  .004.  Other 
parameter  choices  can  result  in  still  faster  small  signal  be¬ 
havior.  A  numerical  simulation  of  a  start-up  transient  using 
these  parameters  is  shown  in  Figure  6.  Note  that  the  initial 
condition  for  the  estimate  of  the  nominal  inductor  was  taken 
as  zero.  • 


Figure  6:  Start-Up  Transient  in  Second  Order  Converter  Us- 
Exan.ple:  Estimation  of  Nominal  Inductor  Current  in  ■  Ad  tive  C<)ntrol  Scheme 
Simple  Up-Down  Converter  In  this  example,  we  apply 
the  adaptive  control  scheme  to  the  secoud  order  up-down  con¬ 
verter  whose  parameters  are  given  in  Section  3.1.  We  now  as¬ 
sume,  however,  that  the  load  is  unknown  but  constant  in  the  More  cmpfe,  examples  that  require  estimation  of  more 
steady  state.  As  a  consequence,  the  nominal  inductor  current  thaa  one  nr,m;n,|  gUte  are  considered  in  [23]. 
is  also  unknown.  This  is  the  parameter  that  our  self- tuning 
mechanism  will  estimate.  In  this  example,  it  is  assumed  that 
the  input  voltage  V,  is  known  (i.e.  measured),  the  nominal 
output  voltage  v„  is  defined  by  the  regulation  problem,  and 

the  nominal  duty  ratio  rf„  is  known  (The  no«mnal  duty  ratio  4  g  ummarizing  Remarks 
can  usually  be  determined  from  V,  and  t>„.)  We  work  with  0 

the  augmented  model 

The  Lyapunov-based  control  described  in  this  paper  is  evi- 
i'  0  (l—dn)/L  0  i  dently  a  promising  approach  to  the  control  of  switched  mode 

(/  =  -(1  -  dn)/C  0  0  v  +  power  converters.  The  method  can  be  extended  to  couvert- 

(i,'„)'  0  0  0  (it*)  ers  that  handle  time- varying  input-output  waveforms,  see  the 

discussion  in  [23].  In  applications  in  distributed  power  sup¬ 
ply  environments,  this  type  of  control  may  prove  very  useful 
since  it  may  become  necessary  to  stabilize  arbitrary  intercon¬ 
nections  of  converters  and  loads.  The  method  also  lends  itself 
y  =  (lj  -  v,)[i -(£«„)} -H,(v  -  u„).  (23)  to  the  design  of  state  observers  as  outlined  in  [23]. 


(V.  -  v,)/L 
«,/C 
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A  Passivity,  Incremental  Passivity, 
and  Relative  Passivity 

In  order  to  state  the  following  definitions  in  a  relatively  gen¬ 
eral  way.  we  assume  the  input-output  vector  pair  of  an  n-port 
to  be  a  hybrid  pair.  That  is,  the  input  u(t)  and  output  y(t)  of 
an  n-port  are  n-component  vectors  whose  elements  represent 
port  voltages  or  currents.  The  components  of  y(<)  are  com¬ 
plementary  to  those  of  u(f),  and  oriented  such  that  u(t)‘y(t) 
is  the  instantaneous  power  entering  the  network  at  its  ports. 
In  the  following,  the  networks  of  interest  axe  assumed  to  be 
time-invariant  unless  otherwise  noted. 

Passivity  The  definition  of  passivity  presented  in  Wyatt 
et.  al.  [17],  Wyatt  [19],  and  in  Hasler  and  Neirynck  [16]  will 
be  adopted  here. 

Definition  A.l  (Available  Energy)  Given  an  n-port  N, 
let  the  available  energy  Eaj  in  state  x  be  the  maximum 
energy  that  can  be  extracted  from  N  when  its  initial  state  is 
x.  with  the  convention  that  E*a  =  +00  if  the  available  energy 
is  unbounded.  That  is, 

Ea*=  Sup/  —n(t)'y(t)  dt  (24) 

T  Jo 

Definition  A. 2  (Passivity)  N  is  passive  if  E^  is  finite 
for  each  initial  state  x. 

Note  that  this  definition  of  passivity  is  directly  tied  to  a  state- 
space  realization  for  the  n-port  in  question.  This  is  not  ob¬ 
jectionable  for  our  purposes  since  we  aim  to  draw  conclusions 
for  switching  converters  for  which  state  models  are  readily 
obtained.  In  the  context  of  a  switching  converter,  the  con¬ 
cept  of  passivity  is  of  use  in  viewing  a  controlled  converter  as 
an  interconnection  of  various  n-ports. 

Incremental  Passivity  The  definition  given  here  follows 
the  svstem  theoretic  framework  of  Desoer  and  Vidyasagar 
[20].  ’ 

Definition  A. 3  (Energy  in  the  Increment)  Given  ann- 
port  N  with  initial  state  x,  let  («i(f),J/i(f))  owi  (uj(f),yi(f)) 
be  any  two  admissible  input-output  trajectories  on  [0,  T\  with 
T  finite.  The  energy  in  the  increment  between  the  two 
trajectories  is  defined  by 

W,{T)  •  [T(«i  -  «a)*(y  1  -  Vs)  dt.  (25) 

Definition  A. 4  (Incremental  Passivity)  An  n-port  N 
with  initial  state  x  is  incrementally  passive  at  state  x  if 
IVt(T),  the  energy  in  the  increment  defined  in  (25),  is  non- 
negatice  for  every  pair  of  admissible  trajectories  on  (0,  T\  with 
T  finite.  If  the  network  is  incrementally  passive  at  all  states 
x  in  the  state  space,  it  is  said  to  be  incrementally  passive. 
The  n-port  is  strictly  incrementally  passive  at  state  x 
if  \Ytl  T)  >  0  whenever  the  two  trajectories  are  distinct.  The 
netwo-i  is  strictly  incrementally  passive  if  it  is  strictly 
incrementally  passive  at  every  state  in  the  state-space. 

Note  tLat  this  definitiou  is  closely  tied  to  the  definition  of  pas¬ 
sivity.  A  passive  network  can  supply  only  finite  energy  while 
an  incrementally  passive  network  can  absorb  only  nonnega¬ 
tive  energy  m  the  increment  between  two  trajectories  (IV,  in 
(25)). 


Relative  Passivity  Incremental  passivity  proves  to  be  too 
strong  a  condition  in  the  case  of  certain  nonlinear  n-ports.  In 
fact,  many  nonlinear  net  works  that  are  not  incrementally  pas¬ 
sive  exhibit  a  closely  related  property  that  we  shall  term  rel¬ 
ative  passivity.  Another  closely  related  notion,  local  passivity 
lot  a  capacitor  (or  inductor)  has  been  introduced  in  [16,  25j. 
However,  our  definition  of  relative  passivity  is  potentially  ap¬ 
plicable  to  any  type  of  network.  To  define  a  relatively  passive 
network,  we  examine  the  energy  in  the  increment  with  respect 
to  a  constant  nominal  operating  point. 

Definition  A. 5  (Relative  Passivity)  Given  an  n -port  N 
with  equilibrium  state  xn  and  nominal  output  y„  correspond¬ 
ing  to  the  constant  input  u,,  consider  the  admissible  trajec¬ 
tory  (u(f),y(f))  on  [0,r]  that  is  obtained  with  initial  state 
x(0)  =  r„.  The  n-port  is  relatively  passive  at  xn  if 

W..(T)  =  jT\u(t)  -  «.]*]y(<)  -  yK)  dt  >  0  (26) 

for  any  finite  T.  The  n-port  is  relatively  passive  if  (26) 
holds  for  any  nominal  operating  point.  N  is  strictly  rela¬ 
tively  passive  at  xK  if  the  inequality  in  (26)  is  strict  when¬ 
ever  x (T)  qt  xm.  N  is  strictly  relatively  passive  if  it  is 
strictly  relatively  passive  for  any  constant  nominal  state. 

In  the  case  of  lossless  elements  for  which  H’r„(r)  is  a  function 
of  only  x„  and  x(T),  H',.( T)  can  be  useful  as  a  Lyapunov 
function. 
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